The kinesin walk: a dynamic model with elastically coupled heads 
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ABSTRACT Recently individual two-headed 
kinesin molecules have been studied in in vitro 
motility assays revealing a number of their pe- 
culiar transport properties. In this paper we pro- 
pose a simple and robust model for the kinesin 
stepping process with elastically coupled Brow- 
nian heads showing all of these properties. The 
analytic and numerical treatment of our model re- 
sults in a very good fit to the experimental data 
and practically has no free parameters. Chang- 
ing the values of the parameters in the restricted 
range allowed by the related experimental esti- 
mates has almost no effect on the shape of the 
curves and results mainly in a variation of the 
zero load velocity which can be directly fitted to 
the measured data. In addition, the model is con- 
sistent with the measured pathway of the kinesin 
ATPase. 



Kinesin is a motor protein converting the energy of ATP 
hydrolysis into mechanical work while moving large dis- 
tances along microtubule filaments and transporting or- 
ganelles and vesicles inside the cytoplasm of eukaryotic 
cells [0. The wall of a microtubule is made up of tubu- 
lin heterodimers arranged in 13 longitudinal rows called 
protofilaments (Fig. |l|). A tubulin heterodimer is 8 nm 
long and consists of two globular proteins about 4 nm in 
diameter: a- and /3-tubulin. The dimers bind head-to- 
tail giving the polarity to the protofilaments. The mi- 
crotubule has a helical surface lattice which can have 
two possible configurations: the adjacent protofilaments 
can be shifted in the longitudinal direction by about 
either 5 nm (A- type lattice) or 1 nm (B-type lattice). 
Electron micrograph measurements of microtubules dec- 
orated with kinesin head fragments 1^-^ indicate the pre- 
dominance of the B-type lattice and show that kinesin 
heads can bind only to the /3-tubulin (meanwhile weakly 
interacting with the a-tubulin also). Native kinesin is a 
dimeric molecule with two globular (~9x3x3 nm) heads. 
Each one has an ATP and a tubulin binding site. 

Recent experimental studies in in vitro motility assays 
have revealed the following properties of kinesin move- 
ment: 

• kinesin moves unidirectionally parallel to the 
protofilaments towards the "(+)" end of the mi- 
crotubule MB: 



• under an increasing load the speed of the kinesin 
decreases almost linearly (see Fig. ^); 

• under its stall force (about 5 pN) kinesin still con- 
sumes ATP at an elevated rate 

• in the absence of ATP (in rigor state) kinesin binds 
to the microtubule very strongly: it supports forces 
in excess of 10 pN Q; 

• the observed step size in the low speed regime (at 
low ATP or at high force) is about 8 nm ||lO(] ; 

• some backward slippage was clearly observed [|| ; 

• the displacement variance at saturating ATP and at 
low load increases linearly with time, but at (only 
somewhat more than) half of the rate of a single 
Poisson stepping process with 8 nm step size, im- 
plying that one step consists of two sequential sub- 
processes with comparable limiting rates pl| . 

Partially motivated by these remarkable experimental 
advances several interesting thermal ratchet type mod- 
els have very recently been developed for the theoreti- 
cal interpretation of the related transport phenomena. 
These models [|3[j20| except that of Ajdari (Ref. [|ll) 
consider cases in which the internal degree of freedom 
of the molecular motors is not taken into account. The 
purpose of the present work is to define and investigate 
a dynamic model for the kinesin walk which has an in- 
ternal degree of freedom and results in a full agreement 
with the experimental data for the range of its param- 
eters allowed by the known estimates for the lower and 
upper bounds of these parameters. 

There are two basic classes of possible models for the 
stepping of kinesin |l^. The first one is the "Long- 
Stride" model in which the heads are moving along a 
single protofilament. The two heads are displaced from 
each other by 8 nm. During the stepping process the 
back head passes the bound front head advancing 16 nm. 
Then the heads change their roles and a new step may 
take place. However, within this framework during the 
long stride the back head has a good chance to bind to 
a (3 site belonging to a neighboring track which is closer 
than the next /3 site in the same track. Another problem 
is that the low displacement variance observed experi- 
mentally cannot be naturally explained. This model was 
studied in detail by Peskin and Oster [Q, who by intro- 
ducing several reaction rate constants found a reasonable 
agreement with the experimental results. 
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FIG. 1. The structure of the B-type hehcal surface lattice 
of the microtubule. The kinesin heads can bind only to the 
/3-tubulin monomers. 

The other possibility is the family of "Two-Step" mod- 
els. These models naturally explain the low displace- 
ment variance because of the two sequential subprocesses. 
During one cycle one of the two heads takes a 8 nm step 
first then the other head steps. Note that the motion 
of just one head advances the centroid of the molecule 
by only 4 nm the successive steps of the two heads in 
rapid pair results in an effective 8 nm step which can be 
observed in the low speed regime. This kind of motion 
gives large stability to the protofilament tracking. The 
differences in the "Two-Step" models arise from the rela- 
tive positions of the heads. If both heads track the same 
protofilament they are not able to pass each other, so 
we can always distinguish the front and the back head. 
In this case the distance of the heads alternates between 
8 and 16 nm, which seems a bit large for the kinesin 
molecule with its 9 — 10 nm long heads, but is still possi- 
ble. (Alternating between and 8 nm is already impossi- 
ble because only one head can bind to each /3-tubulin site 
at a time). Another quite reasonable situation is when 
the two heads track adjacent protofilaments on a B-type 
lattice. In this case the heads are sitting on adjacent /3- 
tubulins displaced from each other by only about 1 nm in 
the longitudinal direction, thus cither can take the first 
step (they can be identical). 



DYNAMICS OF THE KINESIN WALK 

We present here a one-dimensional kinesin walk model 
which describes the whole family of "Two-Step" models. 
Each of the two Brownian heads can move along its own 
one-dimensional periodic potential with period L = 8 nm 
in an overdamped environment. The two potentials can 



be shifted relatively to each other by an arbitrary dis- 
tance (or even can be the same). These potentials rep- 
resent the interaction with the protofilaments and the 
periods are the tubulin hcterodimers. Each period has a 
deep potential valley corresponding to the binding site 
of the /3-tubulin and the other parts of the potential 
are flat. Each valley has an asymmetric "V" shape (see 
Fig. ^) : the slope in the backward direction (towards the 
"(— )" end of the microtubule) is steeper and 0.5 — 1 nm 
long, while the other slope (towards the "(-I-)" end) is 
1.5 — 2 nm long, showing the polarity of the filaments. 
These ranges are in the order of the Debye length. 
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FIG. 2. Schematic picture of the potential and the sub- 
sequent steps (from top to bottom) of the kinesin molecule 
in the two limiting cases. In the case (a) the hinge is always 
in the centroid of the molecule, and advances 4 — 4 nm dur- 
ing both subprocesses. The two heads share the load force 
equally. In the case (b) the relative position of the hinge to 
the back head is fixed, thus the whole load force acts on this 
head. During the first subprocess the hinge does not move, 
however, during the second it advances 8 nm. All the inter- 
mediate cases are possible between (a) and (b). 

In our model the heads are connected at a hinge, and 
a spring acts between them (Fig. H). At the beginning of 
the mechanochemical cycle both heads are sitting in their 
valleys waiting for an ATP molecule, and the spring is un- 
strained. Any configuration of the model is mathemati- 
cally equivalent to a model in which the two potentials 
are identical, the heads are sitting in the same valley, and 
therefore the rest length of the spring is zero. (Further 
on we will consider only this version of the model.) After 
one of the heads binds an ATP molecule, the hydrolysis 
of this ATP causes a conformational change in this head. 
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more precisely, induces the head to take a 8 nm forward 
step. In the language of the model it means that the 
rest length of the spring changes from zero to 8 nm right 
after the hydrolysis. Then, as the first rate- limiting sub- 
process, the strained spring is trying to stretch pushing 
the head to the next valley. Reaching its new 8 nm rest 
length another conformational change occurs in the head 
as a consequence of the ADP release: the rest length 
of the spring changes back to zero quickly, then, as the 
second rate-limiting subprocess, the spring is trying to 
contract pulling forward the other head. Completing 
the contraction a next cycle can start waiting for a new 
ATP molecule. Hirose et al. ^ provide evidence that the 
kinesin-ADP complex has indeed a different conformation 
near the junction of the heads. (Similarly, a myosin head 
also gets over several conformational changes p3| , p^ .) 
An important cooperative feature of our model is that 
only one ATP hydrolysis can occur during each cycle. 

This picture is consistent with the scheme of Gilbert 
et al. for the pathway of the kinesin ATPase. The 
first rate-limiting subprocess was observed as the disso- 
ciation of the head from the microtubule followed by a 
fast rebinding. They did not report on the second one, 
but it is clear that this subprocess already does not be- 
long to the chemical part of the cycle, and induces a slow 
dissociation and a fast rebinding of the other head. 

We can take the load force into account in a natural 
way. Since in the experiments of Svoboda et al. [p|Jlo| 
a large (^^0.5 fim in diameter) and therefore slow (com- 
pared to the kinesin heads) silica bead was linked to the 
hinge of the kinesin molecule by a relatively weak elastic 
tether, we can apply a constant (independent of time) 
force F to the hinge. We arrive at the same conclusion 
by considering the experiments of Hunt et al. H] , where 
a viscous load acted on the moving microtubule while 
the long tail of the kinesin was fixed. But how can we 
transfer the force to the heads when writing down the 
equations of their motion? There are two possible lim- 
iting cases: if the hinge is always in the centroid of the 
molecule (Fig. ||a) the two heads share the load force 
equally, i.e., F/2 force acts on both heads; if the head 
which wants to step is always free from the load force 
(as in Fig. ||b) the whole force F acts on the actual back 
head. Due to the robustness of the model both (and 
therefore all the intermediate) cases yield practically the 
same force- velocity curves (Fig. ||). 

The motion of the heads are described by the Langevin 
equations 



Til = 



F. 



^load ^ ^ 
load 



K 



ix2 - XI - l{t)) + ^l{t) , 
{Xi - X2 - l{t)) + ^2{t) . 

(1) 



where xi and X2 denote the positions of the heads; F 
is the frictional drag coefficient; V{x) is the periodic 
potential; l{t) is the rest length (which alternates be- 
tween and 8 nm due to the conformational changes) 
and K is the stifi'ness of the spring; ^i(i), ^2(0 are 



Gaussian white noises with the autocorrelation function 
(0 W6(i')) = 2kTrS,,,S{t - t') for i,j = 1, 2. 
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FIG. 3. The force-velocity curves for an individual kinesin 
molecule at saturating (^ 90 /iM) ATP in (a) and (c); and at 
low (~10 /iM) ATP concentration in (b). The experimental 
data (scattered symbols) are from Svoboda and Block [9] in 
(a) and (b) obtained by using optical tweezers; and from Hunt 
et al. [8] in (c) applying viscous load. The open circles cor- 
respond to the simultaneous effect of multiple kinesin motors 
in the authors' interpretation [8,9]. The circles and diamonds 
in (a) mean two different set of the measured data. To get 
the best fit (solid lines) we chose 0.7 nm for the backward 
length of the potential valleys and 1.75 nm for the forward 
length; the depth of the valleys was Q = 20kT. In these plots 
the load force was F/2 on each kinesin head. In plot (b) we 
multiplied the average displacement by the ATP consumption 
rate v{catp) = 10 1/s for the best fit. In plot (c) we used the 
same curve as in (a) but multiplied by a factor 0.6, which can 
be explained with the larger viscosity. 
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RESULTS 

In order to obtain results we can compare with the 
experiments first the values of the input parameters have 
to be specified. The drag coefficient F for a single head 
can be calculated from the Stokes formula yielding F sa 
6 X 10-" kg/s. 

The free energy which can be gained from the ATP 
hydrolysis is about 25kT or 100 x 10^^^ J. During a 
stepping cycle the two conformational changes consume 
the free energy of the ATP, while the rest length of the 
spring changes by 8 nm in both cases. This means 
that 1/2 ■ 25fcT « {K/2){8 nmf, from which we get 
K « 1.5 pN/nm. 

A lower limit for the depth Q of the potential valleys 
can be determined from the fact that the kinesin in rigor 
state supports forces in excess of 10 pN. If we try to pull 
out the two-headed kinesin molecule (with drag coeffi- 
cient 2F) from a 2Q deep potential valley with a force F, 
a low limit for the escape rate is 

g-2Q/kT 

kT2T 

as follows from a more general expression (j^) to be dis- 
cussed later. Assuming that the escape rate is much 
larger than 10"^ 1/s we get 13kT as a lower limit of 
Q. In reality the two heads cannot be handled as one 
larger head, because they are not fixed to each other too 
rigidly, therefore, this is a very low limit. Q ~ 20kT is 
expected to be a better estimate. 

Note that the temperature T plays an important role in 
our model due to the deep potential valley the Brownian 
kinesin heads have to escape from. 

We can assume that at the beginning of a stepping 
cycle both heads are sitting in the same valley of the 
potential. After the ATP hydrolysis, as the first confor- 
mational change, the strained spring is trying to stretch, 
pushing one of the heads to a neighboring valley, which 
is L = 8 nm away. If the load force is small there is a 
large probabihty p^^ — J^l/[Jql + "^ol) ^^^^ front 
head jumps to the forward direction due to the asymme- 
try of the potential valleys (the forward slope is less steep 
than the backward), and there is only a small probability 
Pqj^ = JqIj/{Jql + "^ol) ^^^^ the other head jumps back- 
ward. Jq^ and denote the corresponding jumping 
rates. Increasing the load the probability of the forward 
step decreases while that of the backward step increases. 
The average time ioL = 1/('^ol + "^gl) which is needed 
for this stretching also slightly increases. Completing 
this subprocess the second conformational change occurs: 
the spring is trying to contract. Now for low load force 
the probability pjg = J^o/i'^Lo + "^Zo) ^^^^ back- 
ward head jumps forward to the next valley, where the 
other head is sitting, is close to 1, while the probabil- 
ity Plo ~ "^Lo/ i'^Lo + "^Lo) that the forward head jumps 
backward is very small. Increasing the load the situation 
is similar to the previous case. 



Thus, under low load force the kinesin molecule steps 
8 nm forward during almost each mechanochemical cycle. 
But increasing the load the probability that the molecule 
remains on the same place or even takes a backward 8 nm 
step increases. Reaching the stall load the average dis- 
placement of the kinesin becomes zero. 

Since the potential valleys are very deep (compared 
to the thermal energy kT), the heads spend almost all 
of their time at the bottom, and they are able to jump 
only very seldom. Therefore, we can calculate the above 
mentioned jumping rates (J = J^j^, Jq^, J^g or J^^q) 
from an Arrhenius-like formula. All we need for this 
is the effective potential V°^{x) for the jumping head 
at the bottom of the valley and at the top of the bar- 
rier. The effective potential consists of five parts: the 
V{x) periodic potential for the jumping head; the poten- 
tial Fj°^'^x of the external load force; the spring poten- 
tial {{K/2){x - Xr - l{t))'^)^ ; the {V{xr))^ potential 
for the remaining head; and the potential (^Fl°'^'^Xr)^ 
of the external load force for the remaining head. In 
this approximation . means thermal averaging for 
the position Xr of the remaining head at the bottom of 
its valley. Thus, from the corresponding Fokker-Planck 
equation the jumping rate is 

(2) 

where D — kT/T denotes the diffusion coefficient; 
and Vj^fx the minimal and the maximal values of the 
effective potential; and the integrals should be evaluated 
around the appropriate extremum. Due to the depth of 
the potential valleys if we do not average the position of 
the remaining head but assume that it always sits at the 
bottom, we get similar results. 

From the jumping rates one can calculate the aver- 
age displacement d = {PolpIq - PolPlo)^ = (-^ol-^lo " 
JoLJLo){^/i'^OL + JoLW/iJto + Jlo)]L during a cycle 
and the average time t = toL + tw = ^/{Jql + '^ol) + 
1 / ( J^g -I- J^q) needed for this. Increasing the load the av- 
erage displacement decreases due to the increasing prob- 
abilities of the remaining and the backward slippage; and 
the average time slightly increases as a manifestation of 
the so-called Fenn effect (At stall load it is about 

three times larger than without load.) 

At saturating ATP concentration the only rate-limiting 
factor is the stepping process, therefore dividing the av- 
erage displacement by the average time gives the average 
velocity v = d/t of the kinesin (Fig. ||a and c). 

But at low ATP concentration the rate-limiting factor 
is the diffusion of the ATP to the kinesin heads, i.e., 
the stepping process is much faster than getting an ATP. 
Thus the average velocity is proportional to the average 
displacement during one cycle with the prefactor 

/ \ ^sat ' Catp ^ ^sat ^ , 

1^[Catp) — "T; ~ ~ ' Catp ~ const • Catp ; 

-"■m ~r Catp J»-m 
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which is the rate of the ATP consumption. Catp de- 
notes here the ATP concentration; i/gat is the inverse 
average time of a cycle; and is the mechanochem- 
ical Michaehs-Menten constant. Thus the velocity of the 
kinesin is v = v{cATp)d (Fig. |3|3). 

The difference between the shapes of the force- velocity 
curves at saturating (Fig. ||a and c) and at low (Fig. ||b) 
ATP concentration shows unambiguously the Fenn effect. 

Changing the distribution of the load force F between 
the heads causes the change of the average time (toL and 
tio) needed for the two subprocesses, but their sum, i.e., 
the average duration t of the stepping cycle, and the aver- 
age displacement d remains practically unchanged. Thus 
the force- velocity curves are essentially the same for any 
distribution of the load force. 

We have also solved the Langevin equations ([^) by nu- 
merical integration and our values are in full agreement 
with the results obtained from the above analysis. 

DISCUSSION 

Our model describes the stepping process of the ki- 
nesin molecule with elastically coupled heads. We do 
not use reaction rate constants, but only strain induced 
conformational changes and the underlying asymmetric 
periodic potential. Therefore, we can give the full de- 
scription of the dynamics of the motion. 

From Eq. (7) of Ref. it follows that the displace- 
ment variance of a "Two-Step" process (in which each 
cycle consists of two sequential subprocesses with com- 
parable limiting rates) is r = 1 —ps/2 times smaller than 
that of a single Poisson stepping process taking the same 
overall time. The stepping probability Ps < I denotes 
the probability that a cycle produces a forward step. In 
our model at low load force the stepping probability is 
Ps = PoLPto ^ addition, if the ATP concentra- 

tion is saturated the two subprocesses are the only rate- 
limiting factors with the same rates. This leads to r « 0.6 
as compared to r « 0.52 given with an error in the range 
of 0.1 — 0.2 in Ref. jl^. Thus, our model explains the 
low displacement variance at saturating ATP and at low 
load force. Furthermore, fits the measured force-velocity 
curves extremely well, and is consistent with the recent 
experimental studies of the pathway of the kinesin AT- 
Pase. The parameters of the model are essentially deter- 
mined by experiments or theoretical calculations, there- 
fore, they can only be tuned in a rather restricted range. 
The shapes of the force- velocity curves are essentially in- 
sensitive to the parameters within this range, showing 
the robustness of the model. 

If the relative position of the hinge to the backward 
head is fixed as, e.g., in Fig. ^ then only one power 
stroke occurs during one cycle. The main advantage of 
this version is that the hinge advances the whole 8 nm 
distance at once. (Note, that in this case the first sub- 
process is not so important, thus it can be of a different 



type, for example, similar to the myosin step: detach- 
ment, free swinging, and rebinding to the next side.) In 
any other cases (like in Fig. ||a), the two subprocesses 
mean two power strokes, and the hinge advances the 8 nm 
in two parts. In spite of this, one may observe only 8 nm 
displacements, as the two power strokes occur in rapid 
succession. 

From Eq. (||) it follows that the depth and the degree 
of asymmetry of the valleys in the periodic potential we 
assumed have the strongest effects on the results. On the 
other hand, some further details (such as, for example, 
a shallow secondary valley instead of the flat part) have 
less significant influence on the force- velocity curves. The 
temperature plays an essential role in the actual mech- 
anism of the walk. The two major points where tem- 
perature comes into the picture are the following: i) the 
potential valley has to be considerably deeper than 3kT 
to ensure proper binding, ii) on the other hand, from a 
deep valley the heads may escape only if the spring and 
the temperature act simultaneously since the energy of 
the spring cannot be larger than that of the ATP hydrol- 
ysis. 

To decide which theoretical model describes the kinesin 
stepping process in the best way, further experimental 
studies are needed. One possibility is to see how the 
force- velocity curves continue both at higher forces than 
the stall force and at negative forces. Another possibility 
is to measure the displacement variance at higher loads 
and at lower ATP concentrations. This could distinguish 
among the different cases of our model as well, because 
for different distributions of the load force F the average 
duration of the subprocesses are also different resulting 
in different variances. 

We are aware that C. Doering, M. Magnasco and G. 
Oster have considered a similar approach to the problem 
of the kinesin walk . 

The authors are grateful to J. Prost and A. Ajdari for 
useful discussions and for their kind hospitality during 
the visit of I. D. at ESPCI in Paris. The present research 
was supported by the Hungarian Research Foundation 
Grant No. T4439 and No. F17246, and the contract ERB- 
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